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ABSTRACT 

We compare the spectral properties of 227 Gamma Ray Bursts (GRBs) detected by the Fermi Gamma Ray Burst Monitor (GBM) 
up to February 2010 with those of bursts detected by the CGRO/BATSE instrument. Out of 227 Fermi GRBs, 166 have a measured 
peak energy £°^ k of their vF v spectrum: of these 146 and 20 belong the long and short class, respectively. Fermi long bursts follow 
the correlations defined by BATSE bursts between their E*^ vs fluence and peak flux: as already shown for the latter ones, these 
correlations and their slopes do not originate from instrumental selection effects. Fermi/GBM bursts extend such correlations toward 
lower fluence/peak energy values with respect to BATSE ones whereas no GBM long burst with E°^ k exceeding a few MeV is found, 
despite the possibility of detecting them. Again as for BATSE, ~ 5 % of long and almost all short GRBs detected by Fermi/GBM are 
outliers of the -isotropic equivalent energy ('Amati") correlation while no outlier (neither long nor short) of the E p!:ak -isotropic 
equivalent luminosity ("Yonetoku") correlation is found. Fermi long bursts have similar typical values of E°^ but a harder low energy 
spectral index with respect to all BATSE events, exacerbating the inconsistency with the limiting slopes of the simplest synchrotron 
emission models. Although the short GRBs detected by Fermi are still only a few, we confirm that their E°^ k is greater and the low 
energy spectrum is harder than those of long ones. We discuss the robustness of these results with respect to observational biases 
induced by the differences between the GBM and BATSE instruments. 

Key words. Gamma-ray burst: general - Radiation mechanisms: non-thermal 



1. Introduction 

Our current knowledge of the spectral properties of the prompt 
emission in GRBs mainly relies on the data collected in al- 
most 10 years by the Burst And Transient Source Experiment 
(BATSE) onboard the Compton Gamma-Ray Observatory 

\ (CGRO) satellite. 

BATSE was composed by eight identical modules each 
containing two Nal(Tl) scintillation detectors: a Large Area 
Detector (LAD), optimised for sensitivity and directional re- 
sponse, and a Spectroscopy Detector (SD) optimised for energy 
band and resolution. Thanks to its wide energy coverage (~25 
keV - ~1 MeV) and all-sky viewing, BATSE detected more than 
2700 GRBs. The published spectral catalogs comprise relatively 
small sub-samples of bright G RBs, selected on the basis of the 

■ burst fluence and/or peak flux dPreece et al.ll2000l iKaneko et al.l 
l2006h . 

The analysis of such samples revealed two main results about 
the spectral and energetics properties of these GRBs. 

The first one is that the prompt spectra can be typically well 
described by a curved function showing a peak - in a vF v rep- 
resentation - at a typical energy £°^ k of a few hundreds of 
keV but whose distribution spans nearly three orders of magni- 
tude. Large dispersions characterise also the distributions of the 
low- and high-energy photon indices, whos e characteristic val- 
ues are a ~ -1 and yS 2.3, respectively dPreece et al] l2000: 

IKaneko et alj20 06). 
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The second one is that the rest frame peak energy £p e ak of 
long bursts correlates both with the bolometric is otropic energy 
iiiso emitted during the prompt (Amati et al. 2002) and with the 
bolometric isotropic luminosity L p ,j so estimated at the peak of the 
light curve (I Yonetoku et al.1120041) . Such correlations are clearly 
an intriguing powerful clue about the dominant emission mecha- 
nism during the prompt phase. Furthermore, as by correcting by 
the presumed j et opening angle the dis persion of the correlations 
even reduces (Ghirland a et al.ll200 4a). GRBs become potential 
standard candles dGhirlanda e t al. 2004bJ). 

iNava et ail (l2008h (N08 hereafter) performed the spectral 
analysis on a larger fainter sample of BATSE bursts, selected by 
setting a limiting fluence to F = 10~ 6 erg/cm 2 . They found that 
the spectral parameter £°^ k shows a strong correlation with the 
fluence F (or the peak flux P). These imply that the derived dis- 
tribution of £p^ k is strongly affected by a selection based on cuts 
in fluence (or peak flux). The correlations in the observer frames 
(£' obs , —F and E° hs , —P ) may be just the consequence of the rest 

v peak peak ' J J n 

frame correlations mentioned above. Alternatively, it has been 
claimed that the rest fra me correlations are the result of instru- 
menta l selection effects dBand & Preecd 120051: iNakar & Piranl 

I2005h . 

iGhirlanda etafl d2008l hereafter G08) and N08 examined 
two instrumental selection effects which may affect the relation 
between E^^-F and E^ k -P and found that - although they 
do affect the burst sample properties - they are not responsible 
for the correlati o ns fou nd in the observational planes. Moreover 
IGhirlanda et al.l d2010t) recently showed that the correlation 
£peak~ £iso and £p ea k-£p,iso holds for the time-resolved quanti- 
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ties within individu al Fermi bursts (see also iFirmani et alJl2009l 
for Swift bursts and lKrimm et alJ2009l for Swift-Suzaku GRBs) 
and that this "time-resolved" correlation is similar to that de- 
fined by the time integrated properties of different GRBs. These 
results support the hypothesis that the "Amati" and "Yonetoku" 
correlations do have a physical origin and that the trends seen in 
the "observational planes" are just their outcome. 

The same analys is has been applied to short bursts 
dGhirl anda et al] 120091) . though their lower fiuence limits the 
possibility to determine with reasonable accuracy the properties 
of a statistically significant sample of events. Also in this case 
spectral information mainly arise from BATSE data. 

Evidence of a spectral diversity between long and 
short bursts comes from their different hardness-ratios (HR) 
dKouveliotou et al.l [T993). The larger HR of short bursts might 
be ascribed to a larger £°^ k . However, due to the relation be- 
tween £°]^ k and the bolometric fiuence and peak flux, a direct 
comparison between the E° bs , distributions of the two different 

1 peak 

burst classes must ta ke the different fluenc e/peak flux selection 
criteria into account. iGhirlanda et al.l l2009 (G09 hereafter) per- 
formed a detailed spectral analysis of samples of short and long 
BATSE bursts selected on the same limit on the peak flux. They 
found that the peak energy distributions of the two classes are 
similar, while the most significant difference concerns the low- 
energy power-law indices, with short bursts having typically a 
harder a 0.4. 

The existence of intrinsic spectral-energy correlations for 
short bursts is still an open question, due to the paucity of 
short GRBs with measured redshift. However, the findings by 
G09 show that short events do not follow the Speak— Eiso corre- 
lation defined by long ones, w hile their are consistent with a 
similar Epeak-^p.iso correlation (Amati 20061: lOhno et alj feOOS: 
iKrimm et alJl2009h . 

The Fermi satellite, launched in June 2008, opened a new 
and powerful opportunity to shed light on the origin of the GRB 
prompt emission thanks to its two high energy instruments: the 
Large Area Telescope (LAT) and the Gamma-ray Burst Monitor 
(GBM). LAT detected in 1.5 year very high energy emission 
(> 100 MeV) from 15 GRBs, out of which 12 are discussed in 
IGhiseilini et ai](l2010l) . 

The GBM extends the BATSE energy range at both low and 
high energies. Its twelve Nal detectors assure a good spectral 
resolution between ~8 keV and ~1 MeV and two BGO detec- 
tors extend the instrument sensitivity up to ~40 MeV. While - as 
for BATSE - the Nal detectors guarantee full-sky coverage, the 
smaller geometric area (16 times lower than that of the LADs) 
result in a lower sensitivity. 

Until now (beginning of March 2010), about 280 triggers of 
the GBM are reliably classified as GRBo While a spectral cata- 
logue has not yet been published, preliminary results of the spec- 
tral analysis performed by the GBM team have been distributed 
to the community through the Galactic Coordinates Network 
(GCNs) Circulars. For 166 events the spectral properties (E ^) 
are well constrained. 

Their statistically significant number allows to start a mean- 
ingful comparison between BATSE and GBM results. The fact 
that the two instruments have a different energy range and effec- 
tive area sensitivities could be in turn exploited to infer the prop- 
erties of GRBs prompt emission. This is the aim of this work. 
In other words we try to assess whether the GBM is confirming 
any/which ones of the BATSE results or is revealing something 



new about prompt spectra in the hard X/y-ray energy range and 
the existence of the above mentioned spectra vs energetics cor- 
relations. 

In particular the GBM wider energy window should allow 
to examine two important aspects: on one hand the low en- 
ergy photon index - as first inferred from BATSE data - indicate 
(Preec e et al.lll998l 120021: IGhiseilini et al.ll20QQh an inadequacy 
of a simple synchrotron model to account for the GRB prompt 
emission; on the other hand better constraints on the high energy 
spectral index may help to discriminate between quasi-thermal 
and non-thermal processes and would give clues on the energy 
distribution of the emitting electrons. Moreover, due to its wider 
energy range, the GBM could more easily detect soft and hard 
bursts, i.e. those events whose peak energy is outside the BATSE 
energy range. 

In this paper we study these issues with Fermi/GBM bursts. 
In §2 we present the samples of Fermi bursts (both long and 
short) and their spectral parameters, as reported in the GCNs. In 
§3 the properties of Fermi/GBM bursts are compared with the 
BATSE ones in terms of E°\ , F and P , i.e. in the E oh \-F and 

peak' ' peak 

£obs _p pj anes w e compute the relevant instrumental selection 

peak 1 ' 

effects introduced by the GBM characteristics on these quanti- 
ties, considering separately short and long GRBs. We compare 
the spectral parameters of BATSE and Fermi bursts in §4 and dis- 
cuss in §5 the similarities and differences of the spectral proper- 
ties of GRBs between the two instruments and, for each of them, 
between the population of short and long events. 



2. Samples and spectral models 

The spectral ana l ysis performed by different authors 



dPreece et al l l2000t ISakamoto et all [2005; But ler et al.1 120071: 



Nava et al. 2008, K06) revealed that time-integrated and 
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time-resolved burst spectra are (mostly) best fitted by different 
models, the simplest ones being: i) a single power law model 
(PL), ii) a Band function (BAND), which consists of two 
smoothly connected power laws and iii) a Comptonized model 
(COMP), i.e. a power law with a high energy exponential cutoff. 
In some cases a smoothly broken power law with a flexible 
curvature (SBPL) better fits the observed spectra (K06). 

A simple PL function clearly indicates that no break/peak 
energy is detected within the energy range of the instrument. 
However, it is also statistically the best choice when the signal 
to noise is very low (since this model has the lowest number of 
free parameters), as is t he case for most o f Swift/BAT bursts. 

The BAND model (iBand et al.ll 1993b has four parameters to 
describe the low and high power law behaviours, the spectral 
break and the flux normalisation. Typically, the low energy pho- 
ton index a > -2 [N(E) oc E a ] and the high energy photon index 
(3 < -2 [N(E) oc EP~\, so that a peak in vF v can be defined, with 
^peak ran g m 8 between a few keV and a few MeV. 

When there is no evidence for a high energy photon tail or 
^peak ' s near tne hi^h ener gy boundary of the instrument sensi- 
tivity (and /3 is poorly constrained) a COMP model is preferred, 
due to the lower number of parameters. Also in this case a peak 
energy can be defined when a > -2. In the BAND model the 
spectral curvature is fixed by a, (3 and £°^ k - Sometimes it is 
necessary to introduce a fifth parameter that defines its sharp- 
ness, to account for very sharp breaks up to the limit of a broken 
power law function. Analytical expressions and a more detailed 
description of these different spectral models can be found in 
K06. 
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2.1. Long bursts 

For the BATSE instrument we will refer to two different sam- 
ples of long bursts with published spectral properties: the bright 
bursts sample analysed by K06 and the sample of 100 fainter 
bursts published in N08. 

The K06 sampl^ contains all BATSE bursts with peak flux 
P (50 - 300 keV) > 10 photons ctrT 2 s" 1 or fiuence F (~ 20 - 
2000 keV) > 2 x 10~ 5 ergs ctrT 2 . From the total 350 events 
we select the long ones, i.e. with observed duration Tw > 2 s. 
Since we are interested in comparing the distribution of £°^ k , 
we exclude bursts whose best fit model is a PL or for which the 
best fit is a curved model, but with a < -2 or B > -2. This 
selection results in 280 GRBs: for 104 events the spectrum is 
best fitted by a BAND model, for 65 by a COMP model and for 
1 1 1 by a SBPL model. 

The K06 analy sis enlarges previ ous spectral samples based 
on BATSE data dPreece et al.1 12000ft but still contains only the 
brightest objects. For this reason N08 selected a sample of 100 
fainter bursts, representative of the ~ 1000 BATSE bursts with 
fiuence 10~ 6 < F < 2 x 10~ 5 erg ctrT 2 . Of these 44 bursts are 
best described by the COMP model, 44 by the BAND model, 
and 12, not showing evidence of curvature, by a PL function and 
thus are excluded. 

A spectral catalogue of Fermi bursts has not yet been pub- 
lished. However, results of preliminary spectral analysis per- 
formed by the Fermi/GBM team can be found in the GCN 
archive. We collected all bursts with spectral information up to 
the beginning of March 2010 (227 objects) and selected the 146 
of them with known £°^ k and belonging to the long class, with- 
out any cut in fiuence or in peak flux. According to the results 
reported in the GCN, the COMP model better fits the spectrum 
in 90 cases, while in 56 cases a BAND model is required. While 
clearly the analysis reported in the CGN might be in some cases 
preliminary, we expect this does not affect the statistical results 
of our analysis. 



2.2. Short bursts 

The most comprehensive sample of short BATSE GRBs with 
well defined spectral parame t ers is composed by the 79 events 
analysed by iGhirlanda et al.l d2009l) . selected for having P > 
3 photons ctrT 2 s _1 . In most cases their spectrum is best fit with 
a COMP model. 

For the Fermi/GBM instrument we still refer to the spectral 
results reported in the GCN archive. Among the 227 GBM bursts 
with known £°^ k 20 are short (9 are fitted with a COMP model 
and 1 1 with a BAND model). 



3. E obs . -Fiuence and E obs . -Peak Flux planes: 

peak peak ~ 

comparison between BATSE and GBM bursts 

The distribution of GRBs with and without measured redshift 
in the planes E^ k -F and E°^ k -P has been investigated by 
G08 and N08 with the main purpose of probing the role of in- 
strumental selection effects on the two rest frame correspond- 
ing correlations, i .e. the fipeak-Eiso (the "Amati" correlation - 
Ama ti et alJ 2002h and fipeak-^p.iso (the "Yonetoku" correlation 
- Yonetoku et al. H200l . A correlation between the total fiuence 
and £°g^ k was found by Lloyd, Petrosian & Mallozzi (2000) for a 



sample of BATSE bursts without measured redshifts. Thi s find- 
ing has been recently confirmed bv lSakamoto et al.l d2008) using 
a sample of bursts detected by Swift, BATSE and Hete-II. In 
particular, they note that X-Ray Flashes and X-Ray Rich bursts 
satisfy and extend this correlation to lower fluences. N08 consid- 
ered all events with published spectral information detected by 
different instruments (Swift, BATSE, Hete-II, Konus/Wind and 
BeppoSax) together with the 100 fainter BATSE bursts. In both 
planes long bursts define a correlation, with fainter bursts having 
lower £' obs ,. 

peak 

In order to examine the distribution of Fermi bursts in the 
observational planes £ , ° bs , -F and E ohs , -P and compare it with 

r peak nenk r 



peak 
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that of BATSE events it is mandatory to estimate the biased in- 
duced by instrumental selection effects (see G08). One bias is 
represented by the capability of an instrument to be triggered by 
a burst, i.e. the "trigger threshold" (TT). A second one concerns 
the limit imposed on the photon number in order to significantly 
model the spectrum, defined "spectral threshold" (ST) by G08 
and N08. The TT translates into a minimum peak flux, which 
depends on the burst spectrum and in particular on E°^ k and can 

be described as a curve in the E oh \-P plane. The second re- 
peak 1 

quirement (ST) results into a minimum fiuence which not only 
depends on £'° bs , but also on the burst duration. For this reason 

1 peak 

the ST is represented as a region in the E°£* k —F plane. 

Another relevant point is to assess whether Fermi bursts are 
consistent with following the rest-frame correlations defined by 
BATSE GRBs. In the ob servational planes it is possible to define 
a "region of outliers" dNakar & Piranll2005: G08; N08) for the 
corresponding rest-frame correlations (i.e. the "Amati" and the 
"Yonetoku" for the E obs , —F and E ohs . —P planes respectively): 

peak peak r r J ' 

bursts lying in this region, no matter their redshift, cannot be 
consistent with correlations. 

N08 find that 6% of BATSE long bursts are outliers of the 
£peak-Eiso correlation defined by the sub-sample of bursts with 
known redshift, while no outlier is found for the £p e ak-£p,iso cor- 
relation. 

For short bursts the situation is quite different and has been 
investigated by G09. The peak flux selection criterion of their 
sample clearly determines the shape of the distribution of data 
points at low peak fluxes. This fact hamper the possibility of de- 
termining the existence of a correlation, even though a weak in- 
dication for a trend is present (high fluences events, not affected 
by selection effects, tend to have high E°^ k values). 

Interesting results instead regard outliers: almost all short 
bursts are outliers of the £p e ak-Eiso correlation defined by long 
bursts but they can be consistent with their very same £p ea k-£p,iso 
correlation. These findings are supported by the properties of the 
very few short bursts wit h measured redshift with respect to the 
rest-frame correlations (Amati 2006; G09). 



3.1. Estimate of GBM instrumental selection effects 

Following G08, the TT curves are obtained adapting the results 
of Band (2006) and are shown in the right (top and bottom) pan- 
els of Fig.Q] These are clearly the same for long and short bursts 
as the trigger threshold depends only on the peak flux. 

The ST curves have been calculated from numerical simu- 
lations, as done by G08 for BATSE and Swift/BAI. To perform 
these simulations, the typical background spectrum and the de- 
tector response function must be known. For the bursts detected 
by the GBM both of them depend on several factors (e.g. the 
satellite attitude when a burst occurs), and thus there is no uni- 
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Fig. 1. Ej^-Fruence and Ej^-Peak Flux planes for long (upper panels) and short (bottom panels) bursts. Empty squares represent 
BATSE bursts, filled circles Fermi/GBM bursts and filled triangles indicate events detected by other instruments (from N08). In all 
panels the instrumental limits for BATSE and Fermi/GMB are reported: shaded curved regions in the upper left panel show the ST, 
estimated for bursts lasting 5 and 20 s; solid curves in the bottom left panel represents the ST for short bursts. Solid curves in the 
right panels define the TT, identical for short and long events. Thresholds for BATSE are taken from G08 while those for the GBM 
instrument are derived in this work. The dashed curve in the bottom right panel represents the selection criterion applied by G09 for 
their sample of short bursts, i.e. P > 3 phot/cm 2 /s. 



versal background and response matrix which can be adopted. 
To overcome this problem we use the real backgrounds and re- 
sponses of several GRBs detected by the GBM and average the 
results of the simulation to build average ST curves for the pop- 
ulation of long and short bursts, respectively. 

To be consistent with the spectral analysis of the GBM data 
reported in the GCNs, our simulation performs a joint spec- 
tral analysis of the two most illuminated Nal detectors and 
the most illuminated BGO detector. For each detector we ex- 
tract the background spectrum and consider the response files. 
Simulations were perf ormed for the 10 long bursts published in 
Ghirland a et alJ (1201 Oh for two representative values of the du- 
ration, Ttx) = 5 and 20 s (corresponding in Fig. 1 to the curves 
delimiting the red shaded region at low and high fluence, respec- 



tively). For short bursts we estimated the ST curves for all of the 
20 short bursts of the GBM sample assuming a typical duration 
of the simulated spectra of 0.7 s (red curve in the bottom left 
panel of Fig.Q]). This value, also adopted by G08, is the typical 
duration of short GRBs observed by the GBM. 

For the TT and ST curves of the BATSE instrument we sim- 
ply report the results published in G08 (for long bursts) and in 
G09 (for short bursts). 

3.1.1. Long bursts 

^■peak vs Fluence: as can be seen in the top left panel of Fig.[TJ the 
distribution of GBM bursts in this plane (filled circles) displays 
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some differences compared with that of BATSE bursts (empty 
squares). 

The presence of Fermi bursts with low £°^ k (between ~10 
and ~50 keV), not present in the BATSE sample, is clearly due to 
the wider energy range of the GBM instrument, sensitive down 
to ~8 keV (see the ST curves). In this region, GBM bursts are 
consistent with bursts detected by other instruments (filled trian- 
gles). Also GBM bursts define a correlation which mostly over- 
laps with that defined by BATSE bursts and extends to the lower- 
left part of the E oh \-F plane. 

1 peak 1 

Despite the GBM assures good coverage up to ~30 MeV (vs 
~ 1 MeV of BATSE), the distribution of £°^ k at high energies 
is similar for the two instruments. This could suggest that long 
bursts with very high £°^ k (grater than 1 MeV) do not exist or 

are very few. Indeed, we note that high £°^ k also means GRBs 
with high fluences, which are rarer than GRBs with lower flu- 
ences. In fact, the fluence distribution of GBM bursts (Fig. |2l 
empty histogram) is peaked at F ~ 7.8 x 10~ 6 erg cirT 2 , which 
corresponds to typical E^^lOO keV. 

Note that while BATSE bursts appear concentrated at higher 
E^ k -F , the sample is composed by all bursts with F > 2 x 10~ 5 

erg ctrT 2 and only one hundred of fainter bursts representative 
of ~ 1000 objects with 10~ 6 erg ctrT 2 < F < 2 x 10~ 5 erg ctrT 2 . 
The fluence distributions of these two different samples are re- 
ported in Fig. [2] The insert shows the whole fluence distribution 
for BATSE bursts, built up weighing the fainter class by a factor 
10. We refer to this "weighted sample" as "BATSE Total". The 
comparison of the fluence distribution of the BATSE Total sam- 
ple with the GBM one (solid blue and orange lines in the insert 
of Fig. [2j reveals that they are similar (KS probability P = 0.06, 
see Table [1]). 

Another result shown by the GBM bursts and consistent with 
the conclusion drawn from BATSE is that the ST are not re- 
sponsible for the particular distribution of the data in the plane. 
Although they prevent us from detecting bursts with very low 
fluence, they cannot explain why bursts tend to distribute along 
a correlation. This is well visible for BATSE bursts: events with 
large £°^ k tend to concentrate far from the ST, i.e., at higher 
fluences. The trend shown by the BATSE ST (which requires 
higher limiting fluences when is very high and very low) 
cannot explain this behaviour. The same holds for GBM bursts, 
for which the ST are even less curved and cannot be responsible 
for the observed correlation which has a slope ~ 0.20 + 0.04, 
consistent with that defined by BATSE bursts (N08). 

The sensitivity of the GBM instrument also allows to probe 
the region of outliers at high £°^ k - Consistently with BATSE, 
we found that about the 5% of Fermi/GBM bursts are outliers 
with respect to the "Ama ti" correlation defined by the sample of 
GRBs with redshift (see Ghirlanda et al] |2010l for a recent up- 
date). 

£'° bs , vs Peak Flux: the distribution of GBM bursts with re- 
peak 

spect to BATSE bursts is shown in the upper right panel of Fig. 
[TJ Solid curve s represent t he TT derived for both instruments 
(adapted from IB and 120061) . On average, the GBM instrument 
is a factor 3 less sensitive than BATSE in the common energy 
range. As remarked by N08, the sample of BATSE bursts lies far 
from its TT, suggesting that for this instrument the demand of 
performing a reliable spectral analysis dominates the selection. 
This is not the case for the GBM: the data points lie very near the 
TT curves, suggesting that if a burst is detected there is a good 
chance to recover its spectral shape with a certain precision. This 
makes TT and ST competitive selection effects for Fermi/GBM 



Parameter 




BATSE 




Fermi 


KS 




Bright 


Faint 


Total 






Fluence [erg/cm 2 ] 


3.5e-5 


4.9e-6 


5.7e-6 


7.8e-6 


0.06 


[ keV ] 


251 


141 


158 


144 


0.01 



Table 1. Central values of the E° hs , and fluence distributions re- 
peak 

ported in Fig.|2]and Fig. [3] The last column reports the KS prob- 
ability that the distributions for the complete sample of BATSE 
bursts (BATSE Total, see text) and for the Fermi/GBM bursts are 
drawn from the same parent population. 



bursts. As for BATSE, none of the GBM bursts lies in the region 
of outliers: they can be all consistent with the Yonetoku correla- 
tion. 

The £°^ k distribution is shown in Fig. [3] for both instru- 
ments. The central values of the distributions, obtained by mod- 
elling them with a gaussian function, are reported in TableQ] The 
BATSE bright bursts analysed by K06 have a £°^ k distribution 
peaked around 250 keV, at considerably higher energies with re- 
spect to that of GBM bursts which peaks at 140 keV (KS proba- 
bility P = 4 x 10~ 15 ). The difference strongly reduces when con- 
sidering the BATSE Total sample, including both bright (K06) 
and faint (N08) events. The insert in Fig. [3] shows their £°j^ k 
distribution compared with that of Fermi/GBM bursts. They are 
very similar, both in terms of central value and width of the 
gaussian fit (see Table [TJ, suggesting that the two instruments 
are seeing the same burst population (KS probability P = 0.01). 
The only difference regards the low energy end of the distribu- 
tions. For BATSE bursts it has a quite sharp cutoff at ~ 50 keV 
(green hatched histogram in Fig. [3), according to simulations 
performed by G08 (see Fig. [TJ showing that it is very difficult 
for this instrument to recover £°^ k < 50 keV. Large fluences, un- 
usual for such values of E° h \ , would be required. For Fermi the 

peak' ^ 

situation is different: the GBM is sensitive down to ~ 8 keV and 
this should introduce a cutoff on the £ ,obs , distribution at lower 

peak 

values. A sharp cutoff, however, is not visible for two reasons: 
i) the smoother cut introduced by the TT and ST for the GBM 
instrument compared to that for the BATSE instrument (upper 
panels in Fig.QJ); ii) the Fermi sample is not selected in fluence 
(a fluence limit at 10 -6 ergs cirT 2 would deplete the lowest bins, 
producing a less symmetric distribution). 



3.1.2. Short bursts 

For short bursts (bottom panels in Fig. [TJ the situation is quite 
different in both planes. 

Also for the GBM there is only a weak indication of a cor- 
relation in both planes, consistent with the trends suggested by 
BATSE short bursts. However, the overall behavior is different 
from what happens for long GRBs. 

At high £° b ^ k GBM short events occupy the same region of 
the BATSE ones, and even extend the trend E obs ,-F to E° hs , 

' peak peak 

above the BATSE energy range, revealing that £°^ k larger than 
1 MeV are common in short bursts. 

Furthermore, contrary to what one might expect, they appear 
missing below £°^ k ~ 250 keV. At low fluence this can be ac- 
counted for by considering the ST derived for the GBM instru- 
ment: indeed Fig. [TJ shows that GBM short events lie very near 
to the ST curve (bottom left panel of Fig. [TJ, that prevents the 
estimate of £ , °^ k < 250 keV when F<5x 10~ 7 erg ctrT 2 . 
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Fig. 3. E°°l k distributions for BATSE bright bursts (filled his- 
togram, from K06), BATSE faint bursts (shaded histogram, from 
N08) and Fermi/GBM bursts (empty histogram). The insert com- 
pares the E ^, distribution for the BATSE Total sample (see text) 
with that for Fermi/GBM bursts. 



The higher sensitivity of BATSE instead implies that its ST 
is located at lower fluence. This however does not account for 
the absence of short GBM burst at £ ob * k < 250 keV and higher 

fluence (> 5 x 10~ 7 erg cm -2 ): the BATSE sample shows that 
GRBs with low E°^ k and on the right side of the GBM ST do 
exist. Their absence in the present GBM sample seems to imply 
that they are relatively rare and therefore that at a fluence > 5 x 
10~ 7 erg ctrT 2 most of the events should have £°^ k > 250 keV. In 

other words this would support the presence of a E^ k -F (broad) 
correlation. A larger sample of short GBM bursts may however 
reveal some of them. 

In terms of E° b \ this translates into a distribution peaked 

peak 1 

at higher energies compared both to the BATSE one and to the 
GBM long events. 

In the E^ k -P plane similar conclusions can be drawn: for 
both instruments the TT curves (solid lines) do not affect the 
samples. The BATSE one is clearly limited by the selecting cut 
applied on the peak flux (shaded curve), while we remind that 
GBM bursts have not been selected a priori either in fluence or 
in peak flux and their presence requires their detection and the 
possibility of performing a significant spectral analysis. From 
the fact that the peak flux of Fermi GRBs is significantly above 
their TT, we infer that their selection is dominated by the ST. 

Almost all Fermi short bursts lie in the E° h \ -F region of 

peak ° 

outliers and are inconsistent at more than 3<x with the Amati 
relation defined by the sample of 105 long bursts (including also 
10 Fermi bursts). This confirms what found with short BATSE 
bursts (G09). 

On the contrary, only one burst appears to be an outliers of 
the £peak-£p,iso correlation of long GRBs. The hypothesis that 
short and long bursts follow the same £p ea k-£p,iso correlation is 
supported by the few short events with known redshift. 



4. Long bursts: parameter distributions and model 
biases 

Spectral parameters, such as £°^ k and a, show different typical 
values and width of their distributions depending on the fitting 
spectral model (K06). For example, a (E^ k ) derived from the 
COMP model is on average softer (harder) than that derived for 
bursts whose best fit model is a BAND function; the low energy 
power-law index distribution for the SBPL is softer than that de- 
rived from the COMP and the BAND models. These differences 
can be ascribed both to biases introduced by models with fewer 
parameter or by real intrinsic differences between spectra best 
described by different functional forms. 

In this section we quantify them for the BATSE sample 
and then investigate if similar behaviours are found in the 
Fermi/GBM one. We consider only long bursts, as the sample 
of short GBM bursts is still scarce. 

4.1. E ohs , distribution 

peak 

At first, we consider the sample of bright BATSE bursts (K06) 
according to their best fit model (COMP, BAND or SBPL). 

Fig. [4] shows the distribution of £°^ k for each model. The 
BAND and SBPL models provides very similar central values 
«£p°eak,SBPL> * ^ptk.BANo) * 232 keV; KS probability P = 
0.366), while higher E° h ^ k are derived with the COMP model 
(X £ ptakC0MP> a 311 keV). Consistently with the results by K06, 
the KS probability between the distributions from the COMP 
and BAND models is P = 3 x 10~ 4 , while it is P = 1.5 x 10 4 
between the COMP and SBPL models. 

Two different factors can account for the higher £°^ k of 

spectra fitted with the COMP model. When E obs , is close to the 

r peak 

high energy bound of the instrument /3 can be hardly constrained: 
in these cases a COMP model describes the data equally well, 
but it is statistically preferred as it has a lower number of free 
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Model 


a 


BATSE 

0^,25 


peak 


a 


Fermi/GBM 


aeff,25 


peak 


KS 

a 


(GBM vs BATSE) 

E oh \ 

peak 


COMP 


-1.14 


-1.22 


311 


-0.94 


-0.98 


-1.10 


145 


0.004 


1 x 10-" 


BAND 


-0.93 


-1.04 


232 


-0.76 


-0.80 


-0.89 


143 


0.04 


1.5 x 10~ 5 


SBPL 


-1.25 




232 














BAND + COMP 


-1.00 


-1.12 


264 


-0.88 


-0.92 


-1.05 


144 


0.005 


5 x 10~ 12 


ALL 


-1.11 


-1.16 


251 


-0.88 


-0.92 


-1.05 


144 


6 x 10-" 


4x 10" 15 



Table 2. Central values of the distributions of a, a e g and £°^ k - The last two columns give the KS probability resulting from the 
comparison of the BATSE and GBM bursts distributions of a and E ohs ,. 
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Fig. 4. distributions for different spectral models (bright 

BATSE bursts). The three histograms refers to events whose 
spectrum is best fitted by a SBPL function (filled), COMP model 
(shaded) and BAND one (empty). 



parameters. Moreover, the simulations by K06 show that the 
COMP model tends to systematically overestimate the best fit 
E° hs , with respect to the simulated value and thus underesti- 

peak * 

mates the hardness of the low energy spectrum. These effects 
can account for the different E° h \ and a distributions between 

peak 

the COMP and BAND models found here (Fig. Hand Fig. [7]- 
upper panel- respectively). 

Fig.[5]shows the £°^ k distribution for Fermi/GBM bursts, ac- 
cording to the best fit model (empty and filled histograms for the 
BAND and COMP models, respectively). The two distributions 
appear quite similar (KS probability = 0.553), both in terms of 
central values (see Tabled and widths. Contrary to what found 
with BATSE data, there is no indication of a bias introduced by 
the COMP model. 

The KS test indicates a significant inconsistency between 
the £°^ k distributions of bright BATSE and GBM bursts (last 
column in Table |5J. However, when faint bursts (N08) are in- 
cluded, the BATSE Total and GBM sample distributions are to- 
tally consistent, (see Table[T]and the insert in Fig. [3j, indicating 
that the two instruments are seeing the same population of bursts 
in terms of E^- 150 keV. 



Fig. 5. E , distributions for different spectral models 

° peak 1 

(Fermi/GBM bursts). The shaded (empty) histogram refers to 
bursts whose spectrum is best fitted by a COMP (BAND) model. 

4.2. a distribution 



iPreece et all dl998l) and lPreece et all (120021) derive the low en- 
ergy power-law index of a sample of 156 time integrated and 
5500 time resolved spectra of BATSE bursts and discuss its in- 
consistency with the predicted pr operties of the simple st syn- 
chrotron emission models (see also Ghirlanda et al. 20031). 

However they do not consider the values of a derived directly 
from th e fit but the effectiv e low energy spectral index a e s, de- 
fined bv lPreece et ailll998l (also adopted by K06) as 



a e ff - a - (a + 2) 



£fid 

gobs ' 
peak 



(1) 



which represents the effective slope of the fitted spectral model at 
a fiducial energy Em- For the BATSE instrument they set Em 
25 keV. 

In the case of the SBPL model the introduction of a e g is 
not necessary, since the low energy index retrieved from the fit 
already describes the spectral slope at low energies. 

For the BAND and the COMP models, instead, a represents 
an asymptotic value and is always harder than a e g . The differ- 
ence between the two slopes depends (for a fixed value of Em) 
only on £°^ k - For high £°^ k , a and a e ff are very similar, since 
the asymptotic a well describes the true slope of the data in the 
low energy part of the spectrum. When E°^ k is close to £fid=25 
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Fig. 6. Bright BATSE bursts: low energy power law index vs 
£ peak for different be st fit models (SBPL, BAND and COMP, 
respectively). In the middle and bottom panels, filled circles rep- 
resent a derived from the best fit and empty ones a e s estimated 
using Eq.Q] 



keV, a and or e ff differ significantly, since a e ff is strongly affected 
by the spectral curvature. 

In Fig. [6] we show, for all the three spectral models, a as a 
function of E ohs , (filled circles). For the BAND and COMP mod- 

peak x y 

els, also a c g vs E^ k (empty circles) is plotted. As expected, a e s 
is on average softer than a and the difference is more evident 
for low values of E° hs , . Note that this difference is indeed less 

peak 

pronounced for bursts best fit with a COMP model as this re- 
turns on average larger £°^ k (see §4. 1). In the case of the SBPL 

model there is a clear trend between a and £°^ k , revealing that 
the value of a is strongly affected by the energy where the spec- 
trum curves with respect to the low energy bound of sensitivity 
of the instrument. A trend between a and E°^ k is weakly present 
also for burst best fit by the BAND model, and its significance 
increases when considering a e s in place of a. 



Given these systematic effects introduced by the definition 
of Qf e fF and the arbitrary choice of Em which is related to the 
energy range of the instrument, we rely on the fitted a values for 
a comparison between BATSE and GBM bursts. 

At first we consider the BATSE spectra of K06 for the three 
models (upper panel of Fig . |7j . The corresponding central values 
of the gaussian fits are reported in Table [2] We note that for the 
effects discussed above the a distribution for the COMP model 
is softer than for the BAND one (KS probability P = 2 x 10 -4 ) 
and that the SBPL one provides the softest a distribution (KS 
probability P = 5 x 1(T 9 compared with the BAND and P = 
0.046 with the COMP ones). This is also the case when a e g is 
considered (see a e ff,25 in Tabled. 

Note that the fact that the fitted values of a are harder than 
a e g makes the inconsistency of GR B spectra with the limits 
of synchrotron theory even stronger (iPreece et al.lfl998l [2002; 
iGhisellini et al.ll2000t iGhirlanda et al.ll2003l) . 

The bottom panel of Fig. [7] shows the a distribution for 
Fermi/GBM bursts, fitted with either the BAND or the COMP 
model. The central values of the two distributions are reported 
in Tab. [2] Also for Fermi bursts we find that the COMP model 
provides values of a softer than the BAND one, but in this case 
the difference is less significant (KS probability P = 0.017). 
Interestingly the a distribution of Fermi bursts (for both models) 
is harder than that of BATSE bursts: the KS probability that the 
two belongs to the same parent population of is P = 6 x 1 0~ 9 (last 
row in Table|2]i. This probability increases when excluding from 
the sample of BATSE bursts those fitted with a SBPL model 
(P = 0.005 - next to last row). For comparison, we estimate 
the effective low energy index at E^ = 25 keV or e ff,25 (Table 0. 
Also in this case no matter the model its value is systematically 
harder compared to 0^25 estimated for BATSE bursts. For the 
GBM we also report ar e ff,8, i- e - estimated at the low energy bound 
of the GBM instrument (8 keV). For this instrument the differ- 
ence between a and ar e ff,8 is not statistically significant (P = 0.26 
when considering the whole GBM sample). We conclude that for 
the GBM the definition of a e s is not meaningful, as the value of 
a retrieved from the fit is an accurate estimate of the true slope 
at the low energy bound of the spectrum. 

5. Discussion and Conclusions 

In this work we study the spectral properties of Fermi bursts de- 
tected by the GBM instrument up to February 2010. We collect 
a sample of 227 GRBs (37 of which are classified as short) from 
the GCN circulars. For 146 long and 20 short bursts the spec- 
trum could be fitted by a curved model so that the peak energy 
of the vF v spectrum could be constrained. 

The population of long bursts is large enough to allow a 
statistical comparison with the BATSE results. We consider a 
complete sample of BATSE bursts down to a limiting fluence of 
F ~ 10~ 6 erg cm 2 . Fig.|2]shows that the two samples have very 
similar fluence distributions. For both instruments the £°^ k dis- 
tribution peaks around £°^ k ~ 150 keV (Fig. [3j. Despite of its 
larger energy range, the GBM extends the £°^ k distribution of 
long bursts only at low energies with respect to BATSE. The a 
distribution, instead, reveals some difference: GBM bursts have 
on average a harder low energy photon index ((o-gbm) = -0.9 
and (ffBATSE) = -11, KS probability P = 6 x 10~ 9 ). Also for 
GBM short bursts we can draw some conclusions about their 
spectral properties. Their £°^ k distribution is shifted towards 
higher energies compared both to long bursts from the same in- 
strument and to short bursts seen by BATSE. The lack of low 
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Fig. 8. Schematic view of the distribution of long and short GRBs in the E^ k -F and E^ k -P planes. The horizontal dashed line at 
~ 30 keV represents the lower limit for the GBM instrument: the simulations performed in this work show that £°^ k can be hardly 
determined below this value. For the BATSE instrument this limit corresponds to ~ 50 keV. The upper limit for BATSE is at ~ 1 
MeV, while for the GBM there is no upper limit in this plane. The vertical dashed line (left panel) shows an example of fluence 
selection, while the dashed curve (right panel) refers to the photon flux selection criterion adopted by G09. 



energy £°^ k (below ~ 250 keV) can be accounted for by the 
spectral threshold we derived for the GBM instrument (see Fig. 
Q]), while the larger energy coverage allows the detection of £°j^ k 
up to ~ 4 MeV. Fermi data confirm that short bursts have on av- 
erage a harder a compared to long bursts ((ffGBM.shon) 0.46), 

as already found in the BATSE sample by G09. 

Similarly to what has been done for long and short GRBs 
detected by BATSE (N08; G09) we examine the distribution of 
Fermi bursts in the E^^-Fluence and E^^-Peak Flux planes in 
order to study instrumental selection effects and test their con- 
sistency with the rest frame correlations (i.e. £p ea k - £iso and 
£peak - Lino, respectively) defined by GRBs with measured red- 
shifts. 

Our main results are: 

- Long GRBs detected by Fermi follow the same E^^-F and 
E^-P correlations defined by BATSE GRBs (Fig. Q). We 
computed the instrumental selection effects of Fermi/GSM - 
as already done for BATSE (G08; N08): the trigger threshold 
and the spectral threshold are not responsible for the correla- 
tions defined by long GRBs in both planes (see Fig.[T). The 
Fermi/GBM spectral range down to 8 keV with respect to 
the limit of 30 keV of BATSE, allows to extend the correla- 
tions to lower peak energies/fluences. Instead, despite of the 
higher energy sensitivity of Fermi/GBM (up to 40 MeV) no 
long GRBs with £°^ k larger than a few MeV are detected. 
This can be due to a real absence of bursts with such high 
^■peak or to the fact that they have large fluences, thus being 
too rare to be detected during 1 .5 year of Fermi observations. 
We conclude that long GRBs detected by Fermi confirm 
what found with BATSE bursts, in particular that they follow 
a correlation both in the £' obs , —F and in the E ohs , -P plane. 

peak peak r 

Moreover, the fraction of bursts detected by Fermi which are 



outliers at more than 3<x with respect to the £p e ak-£iso corre- 
lation is ~5%, similarly to what found in the BATSE sample, 
while there are no outliers (at more than 3<x) of the £p ea k-iiso 
correlation among Fermi long GRBs. 
- Short GRBs detected by Fermi populate a different region 
in the E^ k -F plane with respect to long events, the firmer 
having a larger peak energy and lower fluences compared to 
the latter ones. This is consistent with what found by BATSE 
and confirms that short GRBs do not follow the "Amati" cor- 
relation but they obey the "Yonetoku" correlation defined by 
long events. 

The comparison of short and a representative sample of long 
BATSE GRBs (selected with a similar peak flux threshold) led 
G09 to conclude that their main spectral diversity is due to a 
harder low energy spectral index in short bursts while their 
of BATSE is similarly distributed. Fermi bursts provide the the 
opportunity of re-examining this result for the population of 
short and long GRBs detected by the GBM and also compare 
their spectral properties with those of the BATSE ones. We find 
that: 



- E 



peak °f snort Fermi GRBs is larger and a smaller that those 
of long one, indicating that short events are harder, both in 
terms of their peak energy and low energy spectral index. 
A comparison between GBM and BATSE short bursts re- 
veals that they have similar a while the E° h \ of short GBM 

J peak 

bursts is larger than that of short BATSE events (see Fig.Q] 
bottom left panel). This information is allowed by the higher 
energies which can be detected by the GBM. Moreover, the 
different E° bs , distribution of BATSE and GBM short bursts 

peak 

is affected by the lower sensitivity of the GBM instrument, 
which misses short bursts at low fluence (and therefore low 
E°\). 

peak 7 
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tal dashed line in Fig. [8}. Therefore, the Fermi/GBM shows that 
short GRBs have larger E°^ k with respect to long ones, contrary 
to what found with BATSE (G09). 

When comparing the £°^ k distribution of short and long 
bursts, different conclusions can be drawn, according to the se- 
lection criterion of the samples. The left panel in Fig. [8] shows 
that a given cut in fiuence (represented by the vertical dashed 
line) would result in different £' obs , distributions between short 

' peak 

and long bursts, resulting from their different location in the 
E^ k -F plane. The right panel in Fig.[8]illustrates, instead, what 
happens for a selection in photon flux. This translates into a curve 
in energy flux: the dashed curve represents the cut applied by 
G09 to select both short and long bursts, corresponding to a 
photon flux larger than 3 ph cm _2 s _1 . This criterion applied to 
BATSE bursts produces similar £°^ k distributions of long and 
short events, as indeed found by G09. The very same criterion 
applied to Fermi/GBM bursts results, instead, in different distri- 
butions, since short bursts can have very high £°^ k values, not 
detected in the sample of long bursts. 




-2.0 -1.5 -1.0 -0.5 0.0 
Low energy photon index 
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Fig. 7. Comparison of the low energy spectral index distribu- 
tions for the according to the spectral best fit model for: BATSE 
bursts (from K06, upper panel) and for the Fermi/GBM bursts 
(this work, bottom panel). For the latter the solid and hatched 
histograms correspond to the BAND and COMP model, respec- 
tively. 



- Fermi and BATSE long bursts have a similar E° h \ while 

° peak 

Fermi events tend to have a harder low energy spectral in- 
dex (Fig. [7]i. 

Fig. |8] shows a schematic representation of the current in- 
formation about the distribution of short and long bursts in the 



E pelk~ F and E pelk~ P P lanes - with respect to BATSE, the GBM 
reveals that long bursts extend to lower E^ k , consistently with 
what previously found with other instruments (mainly Hete-II 
and Swift). 

Despite of the high energy sensitivity, also the £°^ k distri- 
bution of Fermi long events extends only up to ~ 1 MeV. The 
situation is different for short GRBs whose £°^ k reach up to ~ 4 

MeV in the present sample. These high £°^ k were not detectable 
by BATSE, whose sensitivity drops at ~ 1 MeV (upper horizon- 
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